Astrometric Shifts in the OGLE-1 Microlensing Events 

David M. Goldberg and Przemyslaw R. Wozniak 

Princeton University Observatory, Princeton, NJ 08544-1001 

email: (goldberg, wozniak) @astro. princeton.edu 

S ■ ABSTRACT 

a^ 

We measure the astrometric shifts of the light centroids for the microlensing 
events in the OGLE-1 database. Of the 15 consistently detected events, 7 
Q . had significant shifts which we were able to measure with a fair degree of 

QQ . confidence. Those events with large shifts are also expected to be unresolved 

with a background "blend" , and thus, we suggest that we have identified events 
in the OGLE-1 catalog which are strongly blended. Though we concentrate on 
J:^ . the OGLE-1 database, and use the DoPHOT package in order to perform our 

^ . analysis, we suggest that this shift is a generic effect, and should be observable 

^ . in any crowded field. 
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Subject headings: Gravitational Microlensing 



Introduction 



The search for gravitational microlensing events has been very successful. To date, over 
'k>( I 100 candidate events have been detected, primarily by the MACHO collaboration (Alcock 

H : et al. 1993), but also significantly by the OGLE (Udalski et al. 1992), DUO (Alard 1996a), 

and EROS (Palanque-Delabrouille et al. 1997) collaborations. In order to maximize event 
rates, the searches take place toward crowded fields such as the Galactic bulge and the 
LMC (Alcock et al. 1996). 

While searches toward crowded fields allow for a large number of detections, this 
also makes it highly likely that an observed star will have several contributing unresolved 
sources. Since only one of these sources is lensed at any given time, events will, in general, 
be blended (Di Stefano & Esin 1995; Alard 1997; Wozniak & Paczyhski 1997; Han 1997; 
Goldberg 1998). 

For any event, we define the blending fraction, / as: 

f = -^ = — < 1 (1) 
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where Fg is the flux from tlie lensed star, Fi is tlie flux from eacli of the contributing stars 
in the point spread function (PSF), and Fq is the total flux from the observed object. 

If an event is blended, a measurement of the microlensing magnification, Aobs{t), will 
be a systematic underestimate of the true magnification of the source star, A[t): 

A,Ut) = l-f + A{t)f<A{t) . (2) 

For a single lens microlensing event, A[t) evolves in a straightforward way (Paczynski 
1986; 1996): 

A{t) = ,"^!^'t^ ; u{tr = «L„ + (^^)' , (3) 



u{t)^uity + 4: 



where Umin is the impact parameter of the lensed star with respect to the lens in units 
of the Einstein radius, to is the characteristic time of the event, and tmax is the time of 
maximum magnification. Wozniak & Paczynski (1997) show that the parameters { /, Umin, 
to } form an approximately degenerate set, and since / is not known independently, the 
other parameters cannot be known either. 

It has been suggested by one of us (Goldberg 1998) that in very crowded fields a large 
fraction of gravitational microlensing events will exhibit a characteristic astro metric shift. 
If an observed image is composed of two or more true stars, then the center of the PSF 
will be weighted toward the center of light of contributing stars. For an unblended or very 
weakly blended event (/ ~ 1), the center of light is essentially defined by the position of the 
lensed star. However, for a strongly blended event (/ <^ 1) the lensed star may be almost 
anywhere within the PSF and not significantly affect the center of light. 

As the magnification increases, the center of light will shift toward the lensed star. 
Thus, if we consider a blended image which is centered at the origin at baseline flux, and 
then lens a star within the image at (Axq, Ayo), simple geometrical arguments show that 
the center of light will shift toward the lensed star: 



A A Aobs{t) — 1 A A Aobsyt) — 1 f.^ 

This effect has been observed; Alard (1996b) demonstrated a shift in DUO 7^ 2, a 
binary event which was known to be blended photometrically. However, astrometric shifts 
have not been more widely used, and indeed, there is a certain amount of skepticism that 
they may be significant (Han 1997). 
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Goldberg (1998) suggested that by looking at astrometric shift information from a 
micro lensing event, we may approach the problem of blending in two ways. First, those 
events with larger shifts will tend to have a higher degree of blending than those with small 
shifts. Detailed models of a population can potentially yield a loose shift-blending relation. 
Second, by performing a best fit of equations @) for Axq and Aj/o, we can find the position 
of the true lensed star with respect to the observed PSF and the local field. High resolution 
foUowup observations can then be used to determine the brightness of the lensed star, and 
give a direct measure of the degree of blending. 

In this paper, we shall measure the astrometric shift for candidate microlensing events 
in the OGLE archive. Though we do not attempt to explicitly to calculate it here, any 
quantitative estimate of the optical depth of the galaxy to microlensing must include the 
effect of blending, and the presence or absence of a shift can provide a handle on this. Our 
outline is as follows. In § ||, we discuss our analysis methods, in particular, data reduction, 
photometry, and determination of the position shift. In § |^ we discuss the astrometric shifts 
found for the OGLE archive. Finally, in § ^ we discuss some interesting effects that we 
noted in this study as well as present some suggestions for future study. 



2. Analysis 

2.1. Data 

Events analyzed in this paper are those reported by OGLE collaboration from the first 
stage of the experiment (Udalski et al. 1992) including objects caught on the spot by the 
"early warning system" (Udalski et al. 1994b). Wozniak and Szymahski (1997; hereafter 
W& S) extracted from every frame in OGLE-1 archive a 100 x 100 pix subframe centered 
on objects which were formally variable only in one observing season and constant during 
remaining seasons. This dataset includes OGLE-1 microlensing candidates and is restricted 
to good quality frames (grades A-E, F frames rejected). Depending on the field it covers 3 
or 4 seasons between 1992 and 1995, of approximately 4 months of observing each. The 
single Ford (Loral) 2048 x 2048 CCD with 15/im pixels used by OGLE team results in a 
scale of 0.44 arcsec/pix. A 44 x 44 arcsec subframe of a dense stellar field towards the 
Galactic center contains plenty of stars to allow difference photometry and good registering 
of the images. The images, as provided by W& S were de-biased and flat-fielded by an 
automated data pipe line based on IRAF "ccdred" routines. We refer to Udalski et al. 
(1992) and W&S for details of preliminary reductions, selection criteria and occasional 
exceptions. 
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2.2. Method 

All of the remaining steps, that is profile photometry, cross identification of stars and 
registering of images were done independently. Previous data analysis of microlensing 
events has focused on identifying candidate events and using the best photometry possible 
to determine a light curve (and hence the microlensing parameters). The OGLE team used 
a fixed centroid template to "warm start" DoPHOT (Schechter et al. 1995) in order to 
maximize identification rate of the stars in dense fields (equivalent to the assumption that 
there is no shift !). To find the shift, however, we need to define an accurate "absolute" 
reference coordinate system and measure the position of the amplified star with respect to 
this system. Therefore we use DoPHOT in different mode and identify stars in each frame 
independently. 

The outline of our method is as follows. For each of the 100 x 100 pix image, we create 
an observed catalog, that is the list of stars found by DoPHOT. The brightnesses and 
positions of stars we found in all subframes during unlensed period are used to generate a 
template catalog, a basis of the reference coordinate system. During the lensed period, the 
positions of the lensed star relative to the template are computed. Finally, errors in position 
and brightness were determined, and a best fit of equations (^ and (^) were computed. 



2.3. Locating Stars and Photometry 

We benefit from the analysis by the OGLE group in many ways, not in the least 
because we did not need to calculate a seeing and approximate sky level for each subframe, 
as this had already been done. 

We divide the subframes into two categories: those during the event, and those 
preceding or following the event. Again, we have the benefit of the previously computed 
light curves. For those outside the event period (typically \t — tmax\ ^ 2 to), we select one 
of the "A" quality frames to serve as an initial template. Note that this frame is privileged 
primarily in the sense that all other frames will be rotated into the same coordinates. 

The template catalog is adjusted iteratively. We begin with one frame, and 
add additional frames by determining a four parameter coordinate transform: 
{9i, Sxi, 6yi, Fi/ Fi} , that is rotation, parallel shift and fiux ratio. 

We determine the values of the transformation vectors by comparing the new stars to 
the template in order of decreasing brightness. As each star is added, it is rotated using the 
transformation matrix and compared to the previous template, and if the fit is acceptably 
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good (typically Ar < 2 pixels and differences in ffux less than 20%), we iteratively adjust the 
transformation vector to minimize x^ over all the new stars. Otherwise, we throw out the 
star. This continues through the entire catalog. Typically, we were able to associate ~ 120 
stars per subframe with the stars in the template. Once we have found the transformation 
matrix and assigned the new stars to the template stars, we recompute the position and 
brightness of each template star by taking the reference system resulting from an average 
of all the subframes observed thus far. 

We perform this comparison for all the baseline subframes (typically ~ 100), and in 
the end we have well determined centroid positions and brightnesses for the template. 



2.4. Finding the shift 

For the subframes during the event, we want to be careful that we do not use the source 
star or any of its nearby neighbors to compute the transformation vector, since we expect 
a shift as a function of time. As such, we exclude the central 10 x 10 pixels around the 
lensed star for this calculation, but otherwise compute the transformation vector exactly as 
during the unlensed period. Once the transformation matrix is computed we transform the 
position and brightness of the source star to the template coordinate system. 

An estimate of the errors in the shift for each subframe was determined by comparing 
stars similar (differing by less than 20%) to the source star in brightness as measured in 
that frame to their positions and brightnesses in the template catalog. We typically find 
15 — 30 comparison stars per subframe, and their distribution is approximately Gaussian. 
As an estimate of the uncertainties in position, we take the dispersions around the known 
position of each comparison star. 

In order to estimate the uncertainties in brightness, we used the magnitude uncertainty 
given by DoPHOT, multiplied by a correctional coefficient which takes the crowdedness 
into effect (Udalski et al. 1994a). Typically, these coefficients are ~ 1.2. 

Finally, we compute a best ffi to equation @) and, where applicable, equation (|^) Note 
that we only use those data points which occur during the microlensing event. The unlensed 
frames are only used in computing the baseline positions and brightnesses. Uncertainties in 
Axo and Ayo are estimated by: 
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evaluated around the best fit values. 

The uncertainties in the hght curve parameters are computed by determining a 
covariance matrix for Umin, to and tmax] the associated errors are calculated in the standard 
way: 

-, 1 /9 

(6) 



CTi = 


;(c-i)., 


1/2 


where, 

^^ daidaj 


and ctj is the set of light curve parameters 
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3. Results 

Of the 19 candidate microlensing events reported by the OGLE collaboration, 
subframes were reduced for all but OGLE 7^13, which was not part of the normal search 
field, but rather was a foUowup of a MACHO event toward the galactic bulge (Szymahski et 
al. 1994). OGLE ^2 was on the boundary between the BWC and BW5 fields, and hence, 
we have two sets of subframes which are analyzed independently from one another. 

Of the remaining events, OGLE #'s 4, 16, and 19 were not consistently detected using 
DoPHOT. This bears discussion. Since we did not use a template image to determine 
positions of stars, each frame was run independently, and thus dim source stars were not 
always detected. These stars had I magnitudes of ~ 19.3, 18.5, &19.6, respectively. Likewise, 
the lack of a fixed template allowed the centroid to "move" due to local effects. x One of the 
results of this is that, in general, we effectively had worse photometry than in the official 
OGLE catalog. 

We summarize the results of the 15 observed events, including the double set of 
OGLE 7^2 in table 1. Throughout, all units of position are given in units of pixel lengths. 
Of the 15 events, 7 showed shifts of > 0.5 pixels, not including OGLE #12, which will 
be discussed shortly. Some, such as OGLE #6 had such a large and well-defined shift, 
that is straightforward to pick out the slope even by eye. We also show the plots of the 
observed astrometric shift (in pixels) and the light curves for the events in Figures 1 and 
2, respectively. It should further be noted that the microlensing light curves that we 
determined were generally quite similar to those computed by the OGLE collaboration. 
On primary reason for difference in our parameters may have been the use of a moving 
PSF. Since the integrated brightness is essentially a convolution of PSFs, it should be 
unsurprising that allowing a moving PSF will also allow for variation in brightness. 
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Many of the events have interesting, occasionally unexpected, behavior. First, there 
were 3 events which had obvious outliers: OGLE 7^'s 1, 3, and 14. In each case, there were 
points which were several a from the best fit of equation (|^). Moreover, photometrically, 
they clearly differed from the other points in the light curve. We removed the outliers 
only where there was a poor fit both astrometrically and photometrically. For the sake of 
completeness, we have plotted the outliers as open triangles in both Figures 1 and 2. An 
inspection of the frames containing the outliers reveal that they are primarily caused by a 
bad PSF, and CCD defects. 

OGLE ^ 2 was measured twice, and since the shift was significant, (and indeed, 
consistent to within la errorbars in both the x and y directions), this serves as an 
affirmation that we are, in fact, detecting a shift. Moreover, it suggests that our estimates 
on our uncertainties were not unreasonable. 

OGLE #5 had a significant shift of ~ 0.7 pixels. This is quite encouraging as this event 
had been previously considered to be blended based on photometric information (Alard 
1997). 

OGLE 7^7 is a binary event, which the OGLE collaboration determined to be blended. 
They were able to perform a best fit to the binary lens microlensing parameters (Udalski 
et al. 1994c) using an observation of the event at the caustic crossing. We were unable to 
observe this point, and hence our photometry would be insufficient to uniquely determine 
the blending of the event. We did, however, detect a small astrometric shift of ^ 0.16 pixels. 

Near OGLE #12, we found a nearby (~ 3 pixels) dim star which was only occasionally 
observed. When it was not observed, the best fit of the source star was shifted toward the 
dim star. Thus, we find essentially two populations of images, one in which the companion 
is observed (12a), and one in which the companion is not observed (126). Visual inspection 
of the frames reveals a CCD defect near the source image. Since the CCD is not aligned 
exactly the same with respect to the star field from night to night, some observations have 
this defect overlapping the source (in which case, the source is not observed), others have it 
nearby (in which case, the "companion" star is not observed), and others have it far away 
(in which case both are observed). 

The surprising thing about these two populations is that when the companion is 
unobserved, and hence, when the event is more severely blended, we find that the baseline 
brightness of the source is ~ 10% dimmer than when the companion is observed. We believe 
that this is due to the fact that the centroid of the best fit PSF is shifted off the PSF of 
the source star, and hence, a non-negligible fraction of the light is being thrown out by 
DoPHOT. 



The positions for both of the cases are given in coordinates with the weighted average 
center at the origin. Taking the two populations separately, we find that the center of 12a is 
(0.35 , 0.29), and the center of 126 is at (—0.37 , —0.29). Given that, the actual shifts for 
the two cases respectively are (-0.11 ±0.09 , -0.15 ±0.09) and (0.67 ±0.09 , 0.82 ±0.09). 
Thus, 12a has almost no shift. This, coupled with the fact that we have subtracted out a 
nearby star, suggests that the remaining measurement is essentially unblended. 

Note that the values of Axq and Ayo in table 1 are given in the same coordinate 
system, and hence we expect the 12a and 12b to have identical values of both, as Axq and 
Ayo represent the physical position of the source star. Indeed, in the x coordinate, the 
two agree quite well. However, in the y coordinate, the two populations give inconsistent 
values to several a. This is illustrated in Figure Ic, in which we would expect that the shift 
curves of the two observations will intersect at {Aobs — l)/^o6s = 1 (infinite magnification), 
since observed position ought to be the same, whether or not the nearby companion is also 
observed. 

This event is extremely interesting in that it points out the importance of accurate 
centroiding in photometric measurements. Since the PSF centers were fixed in previous 
studies, this has not been an issue. However, if this technique is to be more widely 
employed, the local environment, and the presence or absence of nearby stars must be 
carefully monitored. 

Our light curves for OGLE #'sl4 & 18 are each extremely sparsely populated and of 
very poor quality. This is primarily due to the dimness (1=19, 18.7, respectively) of the 
stars. Though little confidence should be placed in our computed values for the microlensing 
parameters, we have included the curves for the sake of completeness. 



4. Discussion and Conclusions 

The primary result of this analysis seems to be a positive one. Astrometric shifts 
were found convincingly in over half the candidate microlensing events, some, such as 
OGLE 7^6, unambiguously detectable by eye. This confirms the belief that most events 
are significantly blended and that the majority of the events will show some measurable 
shift. Moreover, though we performed analysis on a particular dataset (OGLE-1), using an 
off the shelf image processing software package (DoPHOT), we suggest that the presence 
of an astrometric shift is a generic one, and that all projects toward the Galactic bulge, 
the LMC, or the SMC will show a significant fraction of events with large shifts. Further 
analysis could be used to provide an estimate of a statistical correction to the optical depth 
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of crowded fields to microlensing. 

However, those events which show a shift may also be analyzed using followup analysis 
on a high resolution ground-based instrument. The detection of a shift gives the position 
of the source star in local coordinates with very high accuracy, and hence it becomes a 
straightforward matter to identify the correct source star in a high resolution image, and 
from that determine the value of /, and hence the other microlensing parameters, and the 
correction to the optical depth. Moreover, the very existence of a shift suggest which events 
are expected to be blended. 

However, in doing this, we encountered several difficulties. First, we note that the 
local effects (eg nearby stars of comparable brightness) pose great difficulties in consistently 
measuring the position and brightness of the lensed star. Though this effect was pronounced 
in OGLE #12, it may have played a smaller role in other events and thus positional 
measurements are somewhat problematical. 

Additionally, both local effects, and a dim source star can cause the star to be 
observed only intermittently. As a result, our uncertainties in Axq and At/o as well as the 
microlensing parameters are not measured with as high a precision as if the star had been 
detected on each subframe. 

Moreover, photon statistics, crowding, and variable seeing, and so on, caused some 
difficulty in photometric measurements. Indeed, this is far from the ideal if one wishes to get 
optimum photometry. It should be noted, though, that the photometry does not artificially 
blend additional stars which were not blended in the initial OGLE measurements. If this 
were the case, we would measure a smaller amplification than in the OGLE catalog, which 
we do not. 

One possible (and encouraging) way around these difficulties in measurement is frame 
subtraction (Alard & Lupton 1997). This method will ultimately provide a direct unshifted 
position of the source and unblended light in every frame. 

In conclusion, we have shown that most events toward the galactic bulge do, in fact, 
have a significant astrometric shift. Moreover, we can conclude a high degree of blending. It 
is hoped that followup observations of this and similar analyses can help to provide greater 
accuracy in the measurement of microlensing optical depths. 

The authors would like to gratefully acknowledge Bohdan Paczyhski for invaluable 
discussions, and the OGLE team, in particular Andrzej Udalski and Michal Szymahski for 
allowing us the use of their data. This research was supported by NSF grant AST-9530478, 
and DMG was supported by an NSF graduate research fellowship. 
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^ OGLE ^'s Im 3 & 14 each have 1, 3, and 3 obvious outUers, respectively. In both cases, the outhers 
are both photometrically and astronietrically distinct from the rest of the population. 

^ OGLE #7 is a binary event. 

^ Measurements of OGLE ^12 are highly dependent upon detection of a nearby faint star. The 
interpretation of "events" 12a and 126 are discussed in the text 
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Table 1. A summary of the astrometric shifts and microlensing parameters of the observed OGLE 

events. Column 1 shows the event # used in the OGLE hterature. Column 3 shows the number of 

frames in which the star is observed to be lensed. Columns 4 and 5 show the computed astrometric 

shift (equations 0]) in pixels, and the remaining columns give the computed microlensing 

parameters. Column 8 is the estimated date of maximum amplification in ( J. D. -2448000). 



# 


Star # 


obs. 


Axo 


Ayo 


Araax 
^obs 


to (days) 


trnax (days) 


Xmin 


ll 


BW7 117281 


38 


-0.58 ±0.04 


-0.11 ±0.04 


2.60 ±0.10 


21.3±1.2 


1154.6 ±0.8 


17.6 


2a 


BW5 178651 


26 


-0.39 ±0.04 


-0.53 ±0.04 


8.51 ±0.37 


52.5 ±1.8 


803.3 ±0.6 


21.9 


26 


BWC 10648 


27 


-0.49 ±0.04 


-0.49 ±0.04 


7.53 ±0.35 


48.2 ±1.4 


805.0 ±0.5 


34.7 


31 


BW3 161225 


44 


0.25 ±0.07 


-0.18 ±0.07 


1.30 ±0.02 


14.4 ±1.5 


833.7 ±1.1 


92.8 


5 


BWC 120698 


20 


0.64 ±0.02 


-0.34 ±0.02 


13.58 ±0.16 


12. 7± 0.16 


824.4 ±0.2 


85.1 


6 


MM5-B 128727 


30 


0.86 ±0.06 


0.09 ±0.05 


7.78 ±0.3 


9.3 ±0.2 


818. 7± 0.1 


53.6 


72 


BW8 198503 


29 


-0.04 ±0.01 


0.16 ±0.01 










8 


BW9 138910 


13 


-0.83 ±0.09 


0.93 ±0.12 


1.56 ±0.05 


37.0 ±5.9 


1217.7 ±5.1 


2.4 


9 


MM7-A 86776 


11 


-0.30 ±0.14 


0.76 ±0.15 


1.75 ±0.14 


16.4 ±4.9 


815.0 ±1.5 


6.5 


10 


BW3 161220 


44 


-0.08 ±0.08 


-0.05 ±0.09 


1.11 ±0.01 


103.0 ±21.4 


840.1 ±23.0 


113.1 


11 


BW6 167045 


33 


-0.08 ±0.10 


-0.46 ±0.15 


1.28 ±0.03 


10.4 ±2.0 


1536.8 ±1.9 


16.4 


12a3 


BW5 83758 


27 


0.23 ±0.09 


0.14 ±0.09 


1.94 ±0.09 


17.8 ±1.4 


1583.1 ±1.0 


30.9 


1263 


BW5 83758 


34 


0.30 ±0.09 


0.55 ±0.09 


2.13 ±0.10 


22.1 ±1.4 


1581.6 ±0.9 


74.0 


141 


MMl-A 123474 


5 


-0.5 ±0.15 


-0.43 ±0.15 


2.30 ±0.36 


25.0 ±10.7 


1815.9 ±7.8 


8.4 


15 


BW3 142477 


19 


0.46 ±0.05 


-0.64 ±0.05 


4.81 ±1.68 


16.9 ±0.9 


1854.2 ±0.8 


16.4 


17 


BWIO 176006 


30 


0.00 ±0.10 


-0.05 ±0.07 


1.57±0.15 


139.9 ±52.5 


1975.7 ±55.9 


41.2 


18 


BWl 67895 


6 


-0.07 ±0.25 


3.35 ±0.35 


2.09 ±6.6 


21.2 ±21.7 


1914.5 ±43.3 


0.43 
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Fig. la. — The astrometric shifts in OGLE 4^'sl — 3. The dotted hne shows the best fit to 
equations (4), yielding the values of Axq and Ay^ in Table 1. The solid circle shows the 
position and brightness of the source during the unlensed state. 
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Fig. lb. — As above, the astrometric shifts in OGLE #'s5 — 8. 
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Fig. Ic. — As above, the astrometric shift in OGLE #'s9 — 12. 
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Fig. Id. — As above, the astrometric shift in OGLE #'sl4 — 18. 
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Fig. 2a. — The light curves for OGLE #'sl — 8. The sohd hne shows the newly calculated 
parameters, while the dotted line shows the OGLE parameters. The solid squares are our 
measurements of the magnification. OGLE 7 is a binary, for which we have not calculated 
microlensing parameters. 
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Fig. 2b.— As above, the light curve for OGLE #'s9 - 18. Note, in OGLE 12, the sohd 
squares and sohd hne correspond to OGLE 12a, the "unblended" realization discussed in 
the text. The dotted line and open triangles correspond to the case in which the nearby 
companion was not detected. For OGLE ^'sl4 — 18, the OGLE group had not previously 
pubhshed a light curve. 



